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Abstract 

Our investigations on porous Si show that on increase of pressure it undergoes crystaUine phase 
transitions instead of pressure induced amorphization - claimed earlier, and the amorphous 
phase appears only on release of pressure. This amorphous phase, when subjected to higher 
pressures, transforms reversibly to a higher coordinated primitive hexagonal phase showing a 
kind of memory effect which may be the only example of its kind in the elemental solids. First 
principles calculations and thermodynamic arguments help understand these observations. 
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In the crystalline form, depending upon the size of the crystallites, Si undergoes structural 
modifications to higher coordinated denser forms under pressure. " For example cubic diamond 
phase of nano-crystalline Si (~ 50 nm) transforms directly to the primitive hexagonal (ph) form, 
by-passing several structures observed in the bulk.^ For still smaller sized crystallites, such as in 
porous Si (~ 4-5 nm), the initial cubic phase has been claimed to transform to a high density 
amorphous (HDA) form which on release of pressure was shown to convert to a low density 
amorphous (LDA) structure. Some of these observations have also been supported by 
theoretical calculations. "^ Experiments on bulk amorphous silicon have shown that it 
transforms reversibly to metallic HDA at high pressures . 

In this letter we report two new observations: i) cubic nano porous Si (71-Si) undergoes a 
crystalline phase transition to primitive hexagonal ^ (ph) phase under high pressure and the 
amorphous phase arises only when the pressure is reduced from the stability field of ph phase,( 
Hereafter referred to as nano - amorphous ) ii) this nano- amorphous phase transforms reversibly 
to the ph phase under pressure. We also present theoretical calculations based on density 
fiinctional theory and simple nucleation growth model to explain the experimental observations. 
Moreover, this observed amorphous-crystalline reversible inter-conversion is also observed in 
bulk amorphous Si, though the pressure release needs to be fast unlike in the case of nano - 
amorphous Si. 

The experimental investigations presented here were carried out on 71-Si as well as bulk 
amorphous Si. 71-Si was prepared through electrolytic etching. The particle size, determined 
from x-ray diffraction as well as Raman results (quantum confinement model), was ~ 4-5 nm. 
X-ray structure factor (S(q)) of amorphous Si at ambient conditions was found to be consistent 
with tetrahedral coordination. 99.9% pure amorphous silicon samples were obtained from BHEL 
(India) and were characterized by FTIR and Raman spectroscopy. 

For high pressure experiments the sample was loaded in a -150 \im hole, drilled in a pre- 
indented tungsten gasket of a diamond-anvil cell (DAC). X-ray diffraction experiments on 71-Si, 
were carried out at BLIOXU of SpringS synchrotron, using a x-ray beam of- 30 )j,m (diameter) 
and X= 0.3085 A. 4:1 methanol-ethanol mixture was used as a pressure transmitting medium and 
pressure was estimated using ruby fluorescence method. For amorphous Si, in addition to the 
measurements at SPringS, some of the experiments were also performed at XRDl beamline of 
Elettra synchrotron. For these experiments a few particles of gold were also loaded in the DAC 
to determine the pressure in the cell. Raman measurements on both the samples were carried out 
in our laboratory using micro-Raman set-up. The Raman modes were excited with the solid state 
pumped laser of wavelength 532 nm. 

On increase of pressure, Raman results, shown in Fig. 1, confirm vanishing of Raman 
active peak of cubic phase of 71-Si beyond 18 GPa. For example, at 23 GPa we get a flat 
background. This result is very similar to the results published earlier, which had been 
interpreted as pressure induced amorphization, rationalized in terms of pressure induced pseudo- 
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melting due to negative Clapeyron slope . ' ' However, our x-ray diffraction results, given in 
Fig. 2, show that Ji-Si does not become amorphous at least upto ~ 36 GPa. Instead it transforms 
to 8-coordinated ph phase at ~ 20 GPa. At this pressure the most intense [100] and [101] 
diffraction peaks of the ph phase are clearly visible. These observations are similar to those of 
Tolbert et al^ The coexistence of the diamond as well as the ph phase of 71-Si may be attributed 



to the inhomogeneous distribution of local stresses in the sample, aided by pores etc. Non- 
observability of Raman mode beyond 1 8 GPa is also consistent with this (as ph phase does not 
have any optical phonon) and does not imply an amorphous phase. A small difference in the 
pressure of transformation from that of Tolbert et al^ may be understandable in terms of 
different topology and size of nano-particles. We may also note that in MD simulations for bulk 
Si, cubic diamond phase has been shown to transform to the ph phase on abrupt pressure 
increase . It is also known that P-tin phase becomes inaccessible when Si is compressed at low 
temperatures ( < 100 K).^^ These facts suggest the existence of a low barrier path between cubic 
diamond and ph phase. 

Our Raman results (Fig. 1) show that on release of pressure the primitive hexagonal 
phase obtained from tt-SI transforms to a low density amorphous (LDA) phase through a high 
density amorphous phase (HDA), consistent with the earlier studies.^ The broad band centered at 
392 cm" at ~ 15 GPa , can be assigned to the HDA phase. Emergence of HDA phase in the 
vicinity of the stability region of the B-tin phase implies that the amorphous phase in the present 
case is due to kinetic preference. On further release of pressure the low density amorphous 
phase was identifiable at ~ 4.5 GPa, characterized by a broad hump of tetrahedral Si-Si stretching 
vibrations at 484 cm" . When this LDA phase was re-pressurized, it transformed back to the ph 
phase at ~ 18 GPa, implying reversible amorphous to crystalline phase transformation under 
pressure. From Fig. 1 it is evident that even in the second pressure cycle the LDA-ph and reverse 
transformation proceeds through a HDA phase. In fact when pressure was released in the second 
pressure cycle a coexistence of the LDA and HDA phases was observed at ~ 7.5 GPa. 

Our x-ray diffraction (XRD) results are also in agreement with the Raman measurements. 
On release of pressure the XRD pattern shows that a small fraction of the amorphous phase was 
observable even at 12 GPa. On re-pressurizing the amorphous phase (Fig. 2) the diffraction 
peaks of ph phase could be observed at ~ 1 8 GPa. On full release of pressure the sample stays 
fully amorphous (Fig. 2). Hence both XRD and Raman results confirm the pressure induced 
reversibility between the amorphous and ph phases in n-Si. Here we would like to mention that 
no elemental solid is known to display this kind of memory effect between a crystalline and an 
amorphous phase. 

On subjecting the bulk amorphous Si to high pressure, it also transforms to the primitive 
hexagonal phase, preceded by LDA-HDA transformation, similar to nano - amorphous Si. 
However, in this case the crystallization takes place at a lower pressure of ~ 14 GPa. 
Interestingly, bulk amorphous phase exhibits significant deviation during decompression 
process: if the pressure is reduced slowly, it follows same sequence of transformations as the 
bulk crystalline Si i.e. appearance of )9-Sn phase and subsequent R-8 phase on full release. If 
decompressed fast enough (< ~ 60 sees) we get the amorphous (LDA) phase. Re-pressurizing 
this LDA again takes it back to the ph phase following same transformation sequence i.e., LDA- 
HDA-ph. To check whether the observed reversible amorphous-ph transition is an inherent 
property of amorphous Si or it is specific to the LDA phase obtained from ph phase, we 
prepared another amorphous (LDA) Si sample from bulk crystalline Si following pressure- 
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temperature cycle suggested by Imai et al (with fi-Sn phase as a precursor to LDA-Si). Our 
experiments confirm the same reversible nature of amorphous-ph-amorphous phase transition. 



The preferred crystallization of ph from HDA can be understood in terms of 
thermodynamic arguments used earlier by others to understand pressure induced crystallization 
of glass. For quantitative estimates of pressures of crystallization of competing phases, we used 
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the interfacial energy, equation of state of HDA and of crystalline phases obtained from the 
DFT calculations. Using these, it can be shown that HDA should transform to yff-Sn at ~ 12 
GPa and to the ph phase at ~18 GPa. However, the first principle calculations of ref. 13 show 
that LDA-HDA transformation takes place above 15 GPa. Therefore, HDA transforms 
preferentially to ph phase as the crystallization to the B-Sn phase is prohibited due to the fact 
that it would take place at a pressure which is lower than that of HDA formation. 

To gain further insights, particularly in context of kinetics independent transformation of 
ph phase ( obtained from tt-SI) to LDA, several first principles calculations were also carried 
out on Si nano clusters. DFT calculations were used to optimize ionic positions of small silicon 
clusters of approximately -110 atoms and diameter -1.4 nm. Structural relaxations were 
performed within Generalized Gradient Approximation (GGA) using Projector Augmented 
Wave (PAW) method as implemented in Vienna ab-initio simulation package (VASP) 
starting from spherical clusters. An energy cutoff of 400 eV with gamma point sampling was 
used to find the lowest energy configurations of clusters in a 20 A supercell. The amorphous 
structure was generated by heating cubic silicon cluster to 1500 K and quenching it. Molecular 
Dynamics simulation was then carried out for 0.3 ps at 300 K to equilibrate the cluster and was 
then subsequently optimized at K. The resultant cluster had almost the same volume as the 
diamond structure. 

The structures of the different phases of silicon, obtained after structural optimization, are 
shown in Fig. 3. Our results, indicate that for these small clusters, amorphous phase has lower 
energy than even the diamond structure (difference in Energy/atom = 67.9 meV). We found that 
the P-tin structure does not stabilize in small clusters whereas diamond and simple hexagonal 
structures are retained under structural optimization. This could explain the non-observability of 
the P-tin phase in ji-Si. A similar size dependent crystallization has been observed recently in 

-5-7 

silver. 

In conclusion, this is the first observation of a pressure induced reversible amorphous- 
crystalline transformation in an elemental solid and should encourage further experimental as 
well as theoretical studies. 7i-Si does not amorphize on compression and instead it undergoes a 
crystalline transformation and the amorphous phase arises only on decompression. Moreover, 
our studies show that irrespective of the method of preparation, the amorphous Si always 
transforms to the primitive hexagonal phase under compression. 
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Fig. 1 The Raman spectra of 7i-Si at different pressures. At ambient pressure the Raman peak of 
the diamond phase can be seen at 519.54 cm" . 
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Fig. 2 The x-ray diffraction patterns of 7i-Si in the compression and decompression cycles. The 
background of the empty gasketed cell was subtracted from the diffraction patterns. The pattern 
at ambient pressure shows all the diffraction peaks of the cubic diamond phase of silicon. The 
tungsten gasket peaks have been marked with a * and have also been labeled as 'g'. The 
diffraction peaks of the high pressure phase have been indexed with the 'ph' structure. 
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Fig. 3 Clusters of different structures of silicon as simulated by the first principles calculations, 
a) The cubic silicon cluster has the lowest energy among the stable crystal clusters, b) P-tin 
becomes disordered upon relaxation but retains somewhat higher coordination c) primitive 
hexagonal phase forms a more compact cluster although with a higher energy and d) the 
amorphous structure has almost the same volume as the diamond structure but has the lowest 
energy. 



